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ABSTRACT 

Rapid inflows associated with early galaxy formation lead to the accumulation of self-gravitating gas 
in the centers of proto-galaxies. Such gas accumulations are prone to non-axisymmetric instabilities, as 
in the well-known Maclaurin sequence of rotating ellipsoids, which are accompanied by a catastrophic 
loss of angular momentum (J). Self-gravitating gas is also intuitively associated with star formation. 
However, recent simulations of the infall process display highly turbulent continuous flows. We propose 
that J-transfer, which enables the inflow, also suppresses fragmentation. Inefficient J loss by the gas 
leads to decay of turbulence, triggering global instabilities and renewed turbulence driving. Flow 
regulated in this way is stable against fragmentation, whilst staying close to the instability threshold 
for bar formation — thick self-gravitating disks are prone to global instabilities before they become 
unstable locally. On smaller scales, the fraction of gravitationally unstable matter swept up by shocks 
in such a flow is a small and decreasing function of the Mach number. We conclude counterintuitively 
that gas able to cool down to a small fraction of its virial temperature will not fragment as it collapses. 
This provides a venue for supermassive black holes to form via direct infall, without the intermediary 
stage of forming a star cluster. Some black holes could have formed or grown in massive halos at low 
redshifts. Thus the fragmentation is intimately related to J redistribution within the system: it is 
less dependent on the molecular/metal cooling but is conditioned by the ability of the flow to develop 
virial, supersonic turbulence. 

Subject headings: cosmology: dark matter — galaxies: evolution — galaxies: formation — galaxies: 
halos — galaxies: interactions — galaxies: kinematics and dynamics 



1. INTRODUCTION 

The formation of luminous galactic nuclei and super- 
massive black holes is inherently related to the gravita- 
tional collapse of gas and its concurrent release of large 
quantities of angular momentum. Recent numerical sim- 
ulations of this process (Regan & Haehnelt 2009; Levine 
et al. 2008; Wise, Turk & Abel 2008; Englmaier & 
Shlosman 2004) reveal a number of puzzling empirical 
details. First, the collapse appears to be self-similar, over 
as many as 12 decades in radius, with a density profile 
of p oc when averaged over spheres, where r is the 
spherical radius. This is observed despite an extremely 
asymmetric and chaotic appearance locally. Second, it is 
accompanied by fully developed, highly supersonic tur- 
bulence, with near- virial velocities at all radii. And third, 
it does not show any fragmentation, despite being self- 
gravitating and isothermal. The gas temperature is sta- 
ble at a fewxlO^ K (sound speed Cg ~ 3 km s~^), as 
governed by the atomic cooling rate, and the flow Mach 
number is ~ few. Here we propose an explanation for 
these phenomena. 

Angular momentum transfer is of paramount impor- 
tance to the evolution of many astrophysical systems. 
Such systems, which are kept intact by their internal 
gravitational energy, or self-gravity, are found on all 
scales, from planetary and stellar, to galactic and be- 
yond. It is understood that (long-range) gravitational 
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torques rather than (local) viscous torques are respon- 
sible for angular momentum redistribution in these sys- 
tems (e.g., Lynden-Bell & Kalnajs 1972), especially on 
scales where torques generated by ordered magnetic fields 
(Blandford & Payne 1982) are inefficient. 

The amplitude of gravitational torques in a gaseous 
self-gravitating object is determined by the gravitational 
quadrupole moment, i.e., by its departure from axial 
symmetry. When such an object is supported by ro- 
tation against gravitational collapse, its axial symmetry 
is known to be spontaneously broken when the ratio of 
bulk kinetic energy to (absolute) gravitational potential 
energy, T/|l/F|, is larger than some critical value. The 
Maclaurin sequence of rotating fluid ellipsoids, being an 
angular momentum, J, sequence, can serve as a simple 
quantitative example (see, e.g., Chandrasekhar [1969], 
for a review). At some value of J, the axial symmetry 
of an ellipsoid is broken abruptly by the m = 2 mode to 
a barlike shape, due to the lower energy of this config- 
uration. A more robust stability indicator, which is ap- 
plicable to differentially rotating, centrally concentrated 
fluids, is a = (r/|iy|)/(f7/f7,,), where < a < 1/2, Vl is 
the angular velocity and fij = (27rGpA)^/^ is the Jeans 
frequency, A being the shape factor defined in terms of 
the meridional eccentricity (Christodoulou, Shlosman & 
Tohline 1995). The bar instability appears to be univer- 
sal and is triggered in any self-gravitating system, coUi- 
sional or coUisionless, with various mass density profiles 
and sufficient bulk rotation. The subsequent evolution, 
however, depends on the ability of the system to dissipate 
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energy. 

When self-gravitating systems shed angular momen- 
tum, their potential (i.e., self-gravitating) and total en- 
ergies can decrease as well. The relative amount of angu- 
lar momentum loss, compared to potential energy loss, 
is crucial to determining the outcome. If this ratio is 
low enough, there is again a tendency to break axial 
symmetry and increase the quadrupole moment, as in 
the Maclaurin sequence. In non-dissipative systems, of 
course, the total energy is conserved. Such systems can- 
not be driven into a state of gravitational collapse by the 
loss of angular momentum, because the loss of poten- 
tial energy saturates. On the other hand, in dissipative 
systems, energy loss can keep pace with the loss of an- 
gular momentum, forcing them into continuous collapse. 
The latter situation applies in fluid systems, such as gas- 
dominatod disks, and forms the essence of the 'bars-in- 
bars' mechanism of angular momentum transfer in nested 
stellar/fluid or fluid/fluid bars, which is accompanied by 
dynamical inflow (Shlosman, Frank & Begelman 1989; 
Shlosman. Begelman & Prank 1990). Such fluid bars 
can form in maximally rotating neutron stars (e.g., Shi- 
bata, Baumgarte & Shapiro 2000), trigger the formation 
of black holes (Begelman, Volonteri & Rees 2006), or 
dominate galactic proto-disks at high and intermediate 
redshifts (e.g., Heller, Shlosman & Athanassouala 2007; 
Romano-Diaz et al. 2008), and hence provide a imiver- 
sal channel for angular momentum redistribution in self- 
gravitating systems. 

2. TURBULENCE DRIVING 

To ensure continuous collapse, a rotation-dominated 
self-gravitating system must efficiently decrease its J and 
remain dissipative. Stellar systems (i.e., disk galaxies) 
retain the first property, but plainly lack the second. 
What is less clear is whether gasdynamical systems are 
guaranteed to remain dissipative as they collapse. Being 
dominated by self-gravity, they can fragment and these 
fragments can contract, sharply decreasing their interac- 
tion cross-sections and becoming essentially collisionless. 
Star formation is an example of this process. The in- 
ternal dissipation will be dramatically curtailed by the 
star formation and, consequently, the overall collapse of 
the system will be terminated or at least delayed. How- 
ever, the ability of fragments to contract depends on their 
equation of state and should not be taken for granted 
(Paczyhski 1978; Shlosman et al. 1990). The require- 
ment that inflow accompanies the transfer of J in nested 
fluid bars critically depends on the non-fragmentation of 
the flow. It has never been shown that the gasdynamical 
bars-in-bars mechanism can overcome this limitation. 

Feedback from stellar evolution, e.g., from stars formed 
prior to the inflows discussed here, can also affect the 
fragmentation of the flow. Bromm and Loeb (2003) have 
shown that isolated 2a perturbations in the presence of 
metal-free baryons lead to the formation of a single mas- 
sive clump at the center. When H2 cooling is suppressed, 
these clumps enter the isothermal collapse stage without 
further fragmentation. Dijkstra et al. (2008) argue that 
only a small fraction of high-2 halos will be subject to a 
strong ionizing UV continuum (so-called Lyman- Werner 
background) capable of dissociating the H2 molecules 
and suppressing the molecular cooling and, therefore, 
the gas fragmentation. Furthermore, Population III stars 



can enrich the interstellar and intergalactic medium with 
metals that will allow the cooling to proceed below a 
fewxlO'^K even in the absence of (Begelman et al. 
2006; Wise & Abel 2008, and references therein). The 
focal point of this Letter is that molecular and metal 
cooling is not the determining factor in the fragmentation 
of rapid inflows associated with the galaxy formation. 

The first hints that fragmentation may be avoidable 
in bars-in-bars collapse have been provided recently by 
numerical simulations, which reveal that nested bar in- 
stabilities accompany the gravitational collapse, as pre- 
dicted. The flow encompasses two different dynamical 
regimes: (1) a dark matter (DM)-dominated regime, 
from the growing halo virial radius down to a few pc; and 
(2) a gas-dominated regime, within the central few pc. 
The latter regime is characterized by the dynamical de- 
coupling of self-gravitating gas from the DM background 
(Englmaier & Shlosman 2004) . The characteristic radius 
of a few pc applies to DM halos with virial tempera- 
tures somewhat in excess of 10^ K. More massive halos 
will have correspondingly larger gas-dominated central 
regions, and hence the final product of gravitational col- 
lapse can be much more massive. 

It is not surprising that supersonically inflowing gas 
should become turbulent but somewhat surprising that 
the turbulence should maintain a virial level throughout 
the flow, given its high dissipation rate this requires an 
efficient driving mechanism. On larger scales, the turbu- 
lence can be driven by a variety of mechanisms inherent 
in the hierarchical buildup of structure in the Universe, 
such as mergers. On smaller scales, the turbulence can 
be driven by the gravitational energy released during the 
infall. As gas with some seed J goes into gravitational 
collapse, the centrifugal barrier will slow it down at some 
cylindrical radius R. This will happen concurrently with 
the decay of turbulence and will drive T/|ll^| over the 
instability threshold. The instability develops on a dy- 
namical timescale. The growth of the barlike mode in 
the gas will excite shocks and force the gas to lose rota- 
tional support, forcing it into further gravitational col- 
lapse and pumping energy into turbulent motions. Thus, 
we suggest that the persistence of virial, supersonic tur- 
bulence in a self-gravitating inflow is the consequence of a 
self-regulating instability. This cascade will be quenched 
only when the gas becomes optically thick, resulting in 
the sharp rise of thermal energy compared to turbulent 
energy. We quantify the above processes in §3. 

As an example, we assume a time-independent, self- 
similar, self-gravitating, isothermal disklike configuration 
with Keplerian velocity t^x independent of radius. (The 
latter is a natural outcome for a non-rotating flow [Lar- 
son 1969], but it can be extended to inflows with sub- 
stantial rotation.) The surface density scales as R^^ . If 
the inflow speed is roughly independent of radius (as sug- 
gested by the simulations), it does not affect the radial 
momentum equation (since VY^dvY^/AR « 0). One can 
then construct a crude model for the radial structure 
of the inflow. Neglecting thermal pressure, one finds a 
tradeoff between the typical turbulent velocity Vt and 
the rotation speed v^, v"^ = — 2v^. If Vt is small, 
the rotation is close to Keplerian and the ratio T/|VF|, 
which is roughly proportional to (w^/wk)^, exceeds the 
threshold for bar instability. As the radial inflow slows 
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down, the turbulence decays — to the point where T/|Vl^| 
approaches the stability threshold. The decreasing ra- 
tio Vi/v^ will lead to self-similar vertical structure as 
well, with p oc R~^, as seen in the simulations. (We 
note that the self-similar, self-gravitating disk models of 
Toomre [1982] provide a promising platform for refining 
this picture.) This example also holds, with straight- 
forward quantitative modifications, when the gas inflow 
proceeds in the background potential of the DM halo. 

3. DAMPING FRAGMENTATION IN THICK DISKS 

The next aspect of the flow that needs to be addressed 
is fragmentation. If vt <C and the thermal energy 
dominates over the turbulent one, one can show that the 
disk is strongly unstable to breakup into rings (Toomre 
1964). A direct application of the Toomre criterion, 
which is a generalized Jeans instability applied to razor- 
thin, differentially-rotating disks, shows that the disk be- 
comes more stable as Vt increases at the expense of v^. 
We test whether a self-gravitating gaseous disk which is 
subject to the (global) bar instability is stable locally, i.e., 
whether it is subject to fragmentation. Here, we estimate 
the effect of the finite thickness induced by the turbulence 
on the local stability of the disk. In the following, we re- 
fer to 2-D disks as 'thin' disks, and to disks with a finite 
thickness as 'thick' disks. The latter are fully 3-D rotat- 
ing configurations which are supported vertically (along 
the rotation axis) against gravity by turbulent motions, 
while in the equatorial plane, both rotation and turbu- 
lent motions contribute to this support. 

The dispersion relation for a razor-thin, 2-D gaseous 
disk with a surface density E and a turbulent velocity 
(which replaces the sound speed) is given by 

u"^ = + vlk^ - 2T:Gi:\k\, (1) 

where k is the epicyclic frequency, k is the radial wavevec- 

tor and Lo[k) is the wave frequency in the inertial frame. 
The isothermal disks considered here have a fiat rotation 
curve and consequently = 217^, where Q. is the angular 
velocity in the disk. In an unperturbed self-gravitating 
disk, the local Keplerian velocity in the disk is given by 
w|- = ?;| + 2f^, where = ilR. The tradeoff between the 
turbulent "pressure" and the rotational support against 
gravity, for a fixed vk, is crucial for understanding the 
disk dynamics and its stability. 

Thin gaseous disks are stable to fragmentation if the 
Toomre parameter Q satisfies 

'>-^>\ 

These disks are stabilized by rotation on larger scales 
and by the thermal (or turbulent) pressure on smaller 
scales. Sufficiently large Q ensures that these two scales 
overlap and the disk does not fragment. Both thermal 
and turbulent pressures in these disks are assumed to act 
in the equatorial plane only. 

In reality, a disk with finite, isotropic pressure will have 
a finite vertical thickness, h, which is maintained by the 
turbulence: h = {vt/vK)R- The finite disk thickness 
affects the dispersion relation (eq. 1) by decreasing the 
value of the destabilizing gravity term (which appears 
with the negative sign), as follows: 



2 , 2,2 2^GE|fc| 



The correction factor, (1 + \k\h)^^, has been obtained 
by solving the Poisson equation along the rotation (z) 
axis, assuming that the disk has an exponential density 



profile, T,{z) ~ e 



-\z\/h 



Note that eq. (3) reduces to the 



usual dispersion relation for Jeans instability in the limit 
h/R 1. The mass enclosed within a radius R by an 
isothermal, self-gravitating disk is given by M = 27rEi?^, 
so that by using GM/R^ = v^/R, the dispersion relation 
(eq. 3) can be re-written as 

,,2 
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Vj. 



(4) 



1 + ikR)- 

Introducing dimensionless variables, y = kR, x = 
{v,j,/vK)^ and ( = luR/vk, so that {vi/vk)^ = (1 — x)/2, 
we get 

1 — X o y 



C=2x + 



-y 



(5) 



U)^ = + vfk^ - 



l + \k\h' 



(3) 



2 " l + y( 1^)1/2' 

which is the dispersion relation for self-gravitating, 
isothermal gaseous disks with a finite thickness con- 
tributed by turbulent motions. Our goal is to obtain the 
solution of this equation at the threshold for bar insta- 
bility in such disks. The latter is given by Christodoulou 
et al. (1995) as 

where T = {v^/vk)^{v^R/2G) = x{v^R/2G) is the 
kinetic energy of the disk out to radius R, \W\ = 

J^GMdM/R = v^R/G is the potential energy con- 
tained within the same radius, and / is the shape factor. 
In our dimensionless units, r/|VK| — x/2. The shape 
factor is / = 1 for a razor-thin disk and / = (2/3) (1 -|- 
T/|M^|) for thick, differentially rotating spheroidal mass 
configurations. At marginal global stability of the disk, 
a = ttcrit; by substituting acrit into eq. (6), we obtain 

2^crit ~ 0.57. 

Fig. 1 displays the local and global stability regimes 
in the x — y, i.e., {v^/vk)"^ — kR, plane under imposed 
conditions as described below. The dark-grey regions de- 
lineate regimes of local instability, i.e., that are subject 
to fragmentation. In each frame, gas that lies in the left 
part of the diagram has negligible rotation, which cor- 
responds to the usual Jeans instability. The right part 
of each frame shows rotating turbulent gas, correspond- 
ing rather to Toomre's instability. The left frame de- 
picts stability conditions under the artificial assumption 
that the disk is 2-D (razor-thin). As expected, the right 
hand-side, for x > 0.49, is unstable over a wide range 
of wavelengths, because the turbulent pressure and ro- 
tation are not sufficient to prevent fragmentation. The 
second frame corresponds to a case in which we have (ar- 
tificially) allowed the disk to maintain 5% of its thickness 
compared to what it should actually have based on vt- 
The unstable region has moved to a; > 0.66. When the 
full value of disk thickness is accounted for (right frame) , 
the stability boundary lies at 3; > 0.94, i.e., in the regime 
of rapid rotation. 

The light gray area shows the region of global (i.e., 
bar) instability calculated above, x > Xcrit ~ 0.57. We 
observe that gas fiows that possess fully developed turbu- 
lence must move from left to right in Fig. 1, as the turbu- 
lence decays with time. Hence the gas will first reach the 
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Fig. 1. — Stability of thin and thick disks to fragmentation and to bar formation in the x = (va,/vj<^)-^ — y = kR plane. The dark-shaded 
regions show where the disk is unstable to fragmentation, while the vertical line represents the boundary to whose right the disk becomes 
bar unstable. The left frame describes a razor-thin, 2-D disk; the middle frame corresponds to the case where the disk thickness is only 
5% of the actual thickness; while the third frame corresponds to the actual thickness, kh > 1 above the dashed curve in each panel. 



threshold for bar instability while it is still locally stable. 
When the finite thickness of the disk is accounted for in 
the dispersion relation, the global (i.e., bar) instability 
sets in well before the disk is able to fragment. 

What about fragmentation on small (sub-Jeans) scales, 
due to the compression of gas in the isothermal shocks 
that characterize supersonic turbulence? To address this, 
we adopt a widely used model of supersonic turbulence 
(e.g., Padoan 1995; Krumholz & McKee 2005; for re- 
views: Mac Low & Klessen 2004; Elmegreen & Scalo 
2004; McKee & Ostriker 2007). We estimate the fraction 
of the material swept up by the shocks that becomes 
Jeans unstable, by following the formalism of Padoan & 
Nordlund (2002) and Krumholz & McKee (2005). The 
turbulence is characterized by a lognormal probability 
distribution. 
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with the mean Inx ~ — 0.5iTp, where x = p/pa and po 
is the mean density, and the dispersion Cp « [In (1 + 
3XV4)]0-5. 

Two spatial scales characterize such turbulence: the 
Jeans scale, Aj = {nc^/Gp)'^'^ , and the transition scale. 
As, below which the flow becomes subsonic. The lat- 
ter scale relates to the supersonic velocity dispersion, at , 
by (Tt ~ Cs{R/ Xs)^ , which is valid over a large dynamic 
range; with /3 ~ 0.4 — 0.5 (Larson 1981; Krumholtz & 
McKee 2005). In the overdense regions Aj ~ AjoX~°'^, 
where Ajo is the Jeans length at po- To estimate the 
fraction of swept up gas that is Jeans-unstable over one 
free-fall time, /j, we integrate the lognormal distribution 
from Xcrit = Pcrit/Po = (Ajo/As)^ w 67^ 5^1^ to infinity, 
where we have assumed (3 = 0.5 and 70.5 = 0.5vt/vK ^ 1- 
Strikingly, the dependence of /j on Ai assures that 
fj<,2x 10~^ for 7W ^ 3 — such Mach numbers are 
routinely observed in simulations of infiows associated 
with galaxy formation. This differs from the gas evolu- 
tion in galactic disks, where supersonic turbulence de- 
velops as well, but the turbulent velocities are very sub- 
virial (because of lack of sufficient driving) . Further high- 



resolution numerical simulations are needed to shed light 
on the intricacies of shock dissipation in supersonic tur- 
bulence over a wider range of Mach numbers. 

4. DISCUSSION 

To summarize, we have argued that the hypersonic col- 
lapse of self-gravitating gas driven by nested gaseous bars 
is, counterintuitively, not susceptible to fragmentation. 
Fragmentation is suppressed whenever the gas temper- 
ature remains below the virial temperature; moreover, 
the resistance to fragmentation increases with the infall 
Mach number. If the gas has sufficient time to achieve a 
quasi-equilibrium, when it has cooled down or its turbu- 
lence has decayed, it will be subject to fragmentation. If, 
however, a decrease in the level of turbulence leads to re- 
newed non-axisymmetric instability, as we have argued, 
then the infall will recommence. The slowdown of radial 
inflow in regions of relative stability leads to spikes at 
which a substantially exceeds acrit, leading to the cre- 
ation of a distinct bar and a rapid acceleration of infall. 
Thus, the bars-in-bars mechanism seems to be a highly 
intermittent one, as originally envisaged (Shlosman et al. 
1989, 1990). 

If confirmed by carefully designed numerical experi- 
ments, this scenario could have important implications 
for the formation and evolution of the central regions of 
galaxies and, in particular, the formation of supermassive 
black holes (SBHs). It was suggested that the formation 
of SBH seeds by direct collapse would be suppressed in 
lO'^ K halos that were metal-enriched or had significant 
cooling by II2, because the infalling gas would fragment 
and form stars (e.g., Bromm & Loeb 2003; Begelman et 
al. 2006; Wise & Abel 2008). If this is not true, then 
a much larger fraction of 10** K halos could form single 
massive objects. We note that, within this framework, 
DM halos with larger virial temperatures could form 
more massive SBH seeds via direct collapse, if a sufficient 
amount of gas accumulates in their midst. Moreover, 
this mechanism would open the door to the formation 
or rapid growth of SBHs at later epochs in such massive 
halos, for example, following gas-rich mergers. There are 
also implications for the triggering of starbursts and the 
fueling of AGN in galactic nuclei where an SBH already 
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exists. 
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